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STUDY OF STABILITY OF 
LARGE MANEUVERS OF AIRPLANES 
By 

Emile K . Haddad 

LOCK HEED -GEORG! A COMPANY 


SUMMARY 


A predictive method of nonlinear system analysis is used to investigate airplane stability 
and dynamic response during rolling maneuvers. The maneuver roll -rate is not assumed 
to be constant, and the airplane motion is represented by a set of coupled nonlinear 
differential equations. 

In the first portion of the study the general rolling maneuver is kinematically specified by its 
roll -rate variation p(t) . A method for relating the airplane dynamic response to p(t) is 
developed. The method provides analytical expressions for the motion variables in terms of 
the maneuver descriptor p(t) . The accuracy and utilization of the method is illustrated with 
specific values of airplane parameters. A parameterized family of rolling maneuvers is con- 
sidered, for which the method is used to predict specific dynamic response information, such 
as the dependence of the peak angle -of-attack excursion on the maneuver parameters. The 
accuracy of the method is shown to be substantially better than that of previous linear 
analytical methods based on approximating the maneuver roll-rate by a constant. 

In the second portion of the study the stability and motion of the airplane in response to an 
arbitrary actuation of aileron input is considered. Analytical expressions relating motion 
variables to aileron input are obtained. Explicit analytical bounds on the motion variables 
are derived. A stability criterion which guarantees nondivergence of motion in response to 
aileron actuation is presented. 



INTRODUCTION 


The traditional techniques of flight dynamics analysis have depended on a number of 
simplifying assumptions and approximations, applied to the complete equations of motion of 
the aircraft. These equations may be exactly represented by a set of coupled, nonlinear, 
differential equations. The conventional analytical approaches have generally resorted to 
various ways of approximation aimed at reducing the problem to the solution of a set of linear 
differential equations in which the longitudinal and lateral dynamics are uncoupled. The 
adequacy of such assumptions and approximations leading to decoupling and linearization can 
be justified only for restricted ranges of dynamic performance: small deviations from steady 
trimmed straight flight. 

The loss of dynamic -response information incurred by the conventional simplified 
analyses is forcefully exemplified by the well-known rotational coupling problem. Flight 
tests and simulations have indicated that slender aircraft undergoing rapid rolling maneuvers 
may exhibit certain dynamic response characteristics, such as divergence in yaw and pitch or 
excessive excursions in sideslip and angle of attack, which are not predictable from the 
linearized, decoupled equations. Beginning with the work of W. H. Philips (ref. 1), several 
studies (for example, ref. 2,3,4) were made to investigate analytically the characteristics of 
the airplane dynamic response during rolling maneuvers. All of these studies were based on 
assumptions and approximations sufficient to linearize the equations of motion while retaining 
some degree of coupling . Key assumptions for linearization were the requirements of constant 
roll rate and constant forward speed, negligible second-order terms except those involving 
roll rate, small magnitudes of side and vertical velocities, negligible weight components, and 
linearized aerodynamics. These approximations, by eliminating ail the nonlinear product 
terms, reduce the degree of coupling among the component equations, and consequently dis- 
card some of the dynamic response characteristics attributable to such coupling . Further- 
more, the various assumptions made about the dynamic variables are not always representative 
or realistic . A roll maneuver, for instance, may not be executed at constant roll rate, but 
it would typically involve a gradual increase in roll rate. 

Apart from the analytical studies of the rotational coupling problem, one finds investi- 
gations primarily involving extensive simulations on the analog or digital computer (for 
example, ref. 5). In post-design stages, the computer is quite effective for verifying and 
checking specific information about the dynamic response of a given aircraft. However, the 
computer does not lend itself to the formulation of predictive stability and dynamic response 
criteria in concise mathematical form. 

In this study, the airplane dynamic response and stability during rolling maneuvers is 
investigated using nonlinear system analysis techniques. The investigation is based on a new 
approach for nonlinear system analysis developed by E. K. Haddad at the Lockheed -Georgia 
Company (ref. 6, 7, 8). The study develops two analytical methods for predicting dynamic 
response and stability information from the nonlinear equations of motion representing general 
rolling maneuvers for which the roll rate or the aileron input may exhibit any time variation. 
The first method considers the maneuver to be specified by its roll rate time history p(t) . In 
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the second method, the maneuver is specified by the aileron time-history input 6 (t). This 
method leads to a stability criterion which defines the range of aileron inputs for which 
nondivergent motion can be guaranteed. 
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SYMBOLS 


coefficient matrix in state equations of airplane dynamics x = Ax + u 
matrix representing time-invariant portion of matrix A 
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acceleration due to gravity 

elements of the impulse response matrix H 

Impulse response functions for system representing roll 
equation of motion 

Impulse response matrix for coefficnet matrix A Q/ H(t) = e 
body axes moments of inertia 
mass of airplane 

mapping relating a, p, q, and r to roll rate p 
angular rates about body axes 
average value of p 
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P(t) 


q 

Q(t) 



P-P 

time -history bound on p 

dynamic pressure 

time-history bound on q 

time -history bound on r 

parametrized family of rolling maneuvers 


S 
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parameter defining members of family of 
rolling maneuvers. 

wing area 

time 

control input (forcing function) in state equations of 
motion x = Ax + u 


V 

x 

a 

a 

o 

Aa? 

P 

6 

a 

& 

a 

* 

6 

aq 

6 

ar 

6 

aq 


airspeed 

vector representing motion variables x = [ar # (3, q, r 
ang le of attack 

angle of attack at initiation of maneuver 
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angle of side slip 
aileron deflection 
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convolution operator 

A 
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AIRPLANE DYNAMIC RESPONSE DURING AN ARBITRARY 
ROLLING MANEUVER 


The first portion of this study was devoted to the analysis of the dynamic response of the 
airplane during an arbitrarily specified rolling maneuver described by its roll-rate variation 
p(t). The problem being considered is that of relating the airplane stability and dynamic 
response characteristics during a rolling maneuver to the kinematic behavior in roll p(t). 

In the second portion of the study, the airplane dynamic response is directly related to the 
aileron actuation input 6 a (t). One rationale for relating the airplane dynamics to the roll- 
rate p(t) rather than to the aileron input ^j(t) is the fact that during the preliminary design 
effort the specifics of the flight control system may not be known in detail, yet one may be 
interested in the airplane dynamic performance and capabilities during a rolling 
maneuver, which would then be appropriately specified by the kinematic descriptor p(t) 
rather than the control input ^(t). 


Analytical Method for Predicting Airplane Motion During a 
Rolling Maneuver of Arbitrary Roll-rate Variation 

A main result of this study was the development of an analytical method for predicting 
the dynamic response of the aircraft performing a rolling maneuver of arbitrary roll-rate 
variation p(t). 

The maneuver is specified by the roll-rate time-history p(t), and the method provides 
approximate analytical expressions for the response, in which the roll-rate p(t) appears as 
an unspecified parameter function p(t). 

The method is based on the viewpoint that rolling maneuver may be meaningfully des- 
cribed by specifying the variation of roll-rate p(t). From this standpoint, the first four 
equations of motion may be regarded as a mapping between p(t) and x(t) = [af(t), (3(t), q (t) , 
r (t)] : x(t) = M[p(t)] 
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The method provides,for any given aircraft, approximate analytical expressions for the 
nonlinear functional M, in which the response x(t) is expressed explicitly in terms of the 
argument p(t). The method describes a procedure for obtaining successive expressions repre- 
senting progressively better approximations to the exact mapping M. 

Equations of motion for the “basic case" .- The effort in this study has been primarily 
directed towards consideration of the "basic case" involving a rolling maneuver at constant 
forward speed. An adequate representation of the dynamics of this maneuver is given by the 
five-degree-of-freedom equations of motion shown below: 

a - a + Aa , qSC, a = mg 

o ^ La o 


“ = q-pp“5v- c » <* + v- 


mV or V 


coscp 


~ q ■ pp "^v" C La Aa + V- (cos * " ^ 


P= -r + 


pto + ^( C y/ C Yr r + l/ C Yp P ) +0, o p + ■v sln ' f ’ 


q = 


1 - 1 

2 X 


I - I 


pr /c A a + -l.C •a+^C q^ 

r I y ma 2V ma 2V mq y 


• x y , qSb 
r = pq + 


(V + W C nr r + W C n P P ) 


0 ) 


( 2 ) 

(3) 

(4) 



X 

cp= p 



b r 
2V C -tp P 


wjC, r + C„ c 
2V Ir 


■) 


(5) 

(6) 


9 = q coscp - r sincp (7) 

i|f = q sincp + r coscp (8) 

The equations are referred to principal inertia axes. Rudder and elevator are 
assumed fixed during the maneuver. C n§Q and Cy § are assumed to be negligible. 
The aerodynamic coefficients are linearized. The smiplified nonlinear gravity terms 
in equations (1) and (2), and the kinematic relations (6) through (8), are obtained 
under the assumption that pitch and yaw andles (9 , \|r) remain small during the 
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rolling maneuver. Consequent I/, equations (1) - (6) are sufficient to characterize the 
maneuver. Furthermore, equation (6) can be used to express the bank angle cp as 

'P = Jo t P^ dx * 

This expression can be used to eliminate the variable cp from the nonlinear gravity terms of 
equations (1) and (2). The maneuver would then be representable by the five equations (1) - 
(5) with variables a, (3, p, q, and r. During the initial study of the “basic case” the non- 
linear weight effects were ignored. 

Outline of the method .- Letting x = [or, P, q, rU, the first four equations of motion may 
be written as 

x = A(p(t))x + u(p(t)) (9) 

The system in (9) may be viewed as a nonlinear time-varying mapping from p(t) into x(t): 


M 



x(t) = M[p(t)] 


Note that in (9) the argument p(t) appears in the forcing vector u and the coefficient 
matrix A. Thus, doubling the input p(t) does not, in general, result in doubling the output 
x(t), which makes M a nonlinear system. 
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where the matrix A has been separated into a constant matrix A q and a time dependent 
matrix A. A first approximation to x(t) may be obtained by ignoring the second term in 
(12). Thus, let x](t) be the solution to 


*1 “ A o x i + 


03 ) 


The function x] (t) may be conveniently expressed in terms of the input p(t) by using the 
impulse response matrix (Green's function matrix) of the system (13): 


x i(0 ~ H.(p(t)) 


(14) 


Note that xj(t) takes into account the exact time-history in roll p(t) as it appears in 
u(p(t)), and therefore x^(t) is expected to furnish a better approximation to the response 
x(t) than would be obtained from the previous studies based on the approach of replacing 
p(t) by a constant. 

Referring back to the basic equations of motion in (12), and instead of ignoring the term 
A(p(t))x, we now replace x(t) in A(p(t))x by the first approximation X](t). Let X 2 W be the 
solution to 


= A q x 2 + A(p(t))xj(t) + u(p(t)) 

(15) 

= A q x 2 + A(p(t)) H 1 (p(t)) + u(p(t)) 

(16) 

= A q x 2 + Uj(p(t)) 

(17) 


where Uj(p) - A H^(p) + u(p). Again, x.(t) can be expressed in terms of p(t): 


x 2 (t) = H 2 (p(t)) 

The function x 2 (t) is expected to be a "better" approximation to x(t) than Xj(t). This 
procedure can be repeated, to obtain successive mappings x^(t) = H^(p(t)) which represent 
closer approximations to the mapping x(t) = M(p(t)). 

Analytical expressions for the dynamic response. - The method outlined above may be 
used to obtain explicit analytical expressions for the dynamic response successive approxima- 
tions x^(t) In terms of the maneuver roll -rate p(t). For this purpose, let H(t) be the impulse 
response matrix corresponding to the coefficient matrix A 0 : 


H(t) = 


Aot 


(18) 
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Repeating the steps outlined above, and using H(t) to relate x^t) and p(t), one obtains the 
following general expression for x^(t): 

*n W = jj =) / 0 u k (f ' z)dz (l9) 

Where the terms u^(t,z) are defined by the following recurrence relationship 



H(t-z)u(z) 

H(t-z)A(p(z)) f u(z,y)dy 
° k-1 


k= 1 

k = 2, 3, — n 


( 20 ) 


Accuracy of the Method and the Analytical Expressions 

The useful Iness of the above predictive method and the analytical expressions derived 
from it depends on two considerations: 

1 . Does the sequence of approximate expressions jx^(t) j obtained by the method 
converge to the exact airplane response x(t)? 

2. If convergence is guaranteed, how fast is the rate of convergence? What is the 
magnitude of the error expected in the successively closer approximations./ Xj(t), x-^Ct) , Xg(t), 


}• 

Convergence of the analytical expressions | x^(t)| the exact airplane response x(t). - 

The Following is an analytical proof which establishes rigorously the convergence of 
{x n (t)j to x(t): 

Recall from (9) and (12) that the aircraft response during a rolling maneuver p(t) is given by 
the solution of the equation 


x = A(t)x + u(t) 


( 21 ) 


x = A x + A(t)x + u(t) 
o 


( 22 ) 
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and from (19) and (20) the expression for the nth approximate response is 

n t 


x - I f u. (t,z)dz 
n k=l 0 k 


( 23 ) 


where 


H(t-z) u(z) 


k = 1 


u k (t 7 z) 


H(t - z) A(z) j u, ,(z,y)dy k = 2 , 3, . ..n 
At 0 


(24) 

(25) 


H(t) -e 


We now shov/ that is indeed the exact solution to the airplane response equation (22). 

00 \ 

x oo = Z £ u k (t 7 z)dz (26) 


00 t 

= I _d_ J* u, (t,z)d: 

I . A K 


k=l dt 0 

Using Leibnitz rule for differentiation, one has 

oo r t 


. °° i-r 

x oo = j^if A + u k (t # t) 


From the expression for H(t), one obtains 


v 

H= A e 0 = A H 

O 0 


H(0) = I (Identity Matrix) 


From (24) and (25) one has 


Uj(t,t) = H(0) u(t) - u(t) 


~ t 


u k (t,t) - A(t) ^ u k _ 1 (t / y)dy 


k=2,3,. .. 


(27) 


(28) 


(29) 



( 31 ) 


•Jp U k(t/z) = A Q H(t - z) A(z) ^ \_]( z / y /<ty 


= A o u k (t,z) 

Rewriting the expression in (27) as 


k = 2,3, 


oo p- ]■ w 

k=l li » u k (, ' z>dz ] + °1 (m) + k=2 "kM 

and substituting the expressions (28) through (31), one obtains 

00 f GO f 

x = I J A u, (t,z)dz +u(t) + Z A(t) f u (t,z)dz 

“ k=l 0 ° k k=2 0 k-1 


. 00 t 

X = A I 


00 t 


oo o 


J* u k (t,z)dz + u(t) + A(t) I ^ u k _ 1 (t,z)d: 


. 00 t 

x = A I 


00 t 


oo o 


(32) 


L u^(t,z)dz + u(t) + A(t) Z J u^(t,z)dz 

Substituting the expression from (26) into (32), 

• ^ 
x = A x + A(t)x + u(t) 

00 Q 00 00 

Thus x is the exact solution to (22), and therefore represents the required airplane response, 


Accuracy of the approximate responses * n (0 for a parametrized family of rolling 

maneuvers. - A study was conducted in order to assess the comparative accuracy of the 
approximations x], X2, and xg resulting from different members of a class of rolling 
maneuvers. A second objective of the study was to compare the accuracy of xj , x 2 , 
and X3, as obtained by the present method, to the accuracy of the approximation x c (t) 
obtained by replacing p(t) in the equations of motion by a constant (average value of 
p(t) over the duration of the maneuver). The response x c (t) represents the type of result 
that would be obtained if one attempts to apply previous methods (Stemfield, Moul, 
and Brennan) for the prediction of the airplane dynamic response to a realistic rolling 
maneuver comprising a gradually increasing roll-rate followed by a gradually decreasing 
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roll-rate. The study was carried out using the numerical values of the mass and aerodynamic 
parameters of the F-100A. These parameters are shown in Table I . The parametrized class of 
rolling maneuvers used is given by 

-0.6t 

p =P (1 - e ) o < t < 2 tt 

r m - - tt- 

m P 

m 

1 *2TVp A 

= P (e P m-l)i°- 6t t >2 tt 

m -p— 

m 



Each member of the family |p(t)| 
P 


p , corresponding to a specific choice of the parameter 
r m 


m , represents a banking maneuver of 360° executed with a different roll-rate variation. 
The value of P m is indicative of the level of aileron applied during the execution of the 
maneuver. The study was made using the following values of P m : 


P — 1/ -1.5/ - 2 . 0 , -2.5, - 3 . 0 , -3.5, -4.0 rad/sec, 
m 


For these values of P m , the members of |p(t)| are shown in figure 1 . 


It should be noted that, for the configuration parameters used in this study, the Phillips 
criterion indicates the airplane is unstable for roll-rates between -1 .86 and -2.33 rad/sec. 
This implies that, for the maneuvers represented by the values P m = -2.5, - 3 . 0 , -3.5, -4.0, 
the airplane goes thru unstable modes of motion during the execution of the maneuver. This 
is depicted graphically in figure 1 , which shows a plot of |p(t)j.p m . 

For each value of P m , the computer was used to obtain the response x(t), the 
approximate responses x](t), X 2 O), xg(t) as provided by the present method, and the approxi- 
mate response x c (t) based on the constant roll assumption of previous methods. For each 
value of P m , the time histories of a, a], 0 g, 0 ( 3 , a c , 8 , 8 | , $2> £ 3 # and P c were plotted. 

The results obtained are shown in figures 2 thru 15. Inspection of these results leads to the 
following conclusions: 

1 * Although all of the roll maneuvers considered amount to the same bank-angle, 
namely, 360°, the resulting time-histories of a and (3 are strongly dependent on the variation 
of the roll-rate p(t). This fact points to inadequacy of replacing the roll-rate p(t) by a 
constant value for realistic maneuvers involving variation in roll-rate. 
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TABLE I 


MASS AND AERODYNAMIC PARAMETERS USED IN STUDY 


m, slugs 

745 

• x / slug -ft 2 

10,976 

2 

1 , slug -ft 

y 

57,100 

2 

1^, slug-ft 

64,975 

C y , per radian 
3 

-0.28 

Cy , per radian 
r 

0.34 

Cy , per radian 
P 

0.15 

C. , per radian 

a 

3.85 

C , per radian 
nr 
a 

-0.36 

C , per radian 
m 

q 

-3.5 

C , per radian 
m • 

a 

-1 .25 

C , per radian 

0.057 


3 


C , per radian 

-0.095 

r 

C , per radian 

-0.034 

P 

C. , per radian 

-0.255 

P 

C . , per radian 

0.09 

r 

C per radian 

0.06 

C . , per radian 

-0.044 

6 

a 

V, ft/sec 

691 

s, ft 2 

377 

c, ft 

11.3 

b, ft 

36.6 

q, Mt 2 

197 


16 


Ksec) 



2. The peaks of ^occurred consistently during the termination phase of the maneuver. 
The peaks of p occurred during the termination phase of the maneuvers represented by P m = 
-1/ -1 .5. Relatively large peak values of (3 still occurred in the termination phase of the 
maneuvers P m ^ 2. Again, these facts cannot be predicted from an analysis based on a 
constant rol l-rate approximation . 

3. The error plots in figures 16 thru 22 indicate a fairly good rate of convergence of 
the approximate responses x n (t) towards the response x(t) . 

4. The results in figures 16 thru 22 indicate the substantially better results obtained 
by the present method over the methods based on approximating the maneuver roll-rate by 
a constant value. For each of the maneuvers considered, the response x^(t) represented a 
significantly better approximation to the airplane response x(t) than the approximation 
x c (t) obtained by replacing p(t) by its average value. 

Comparison with constant roll-rate methods .- In the previous section, comparison was 
made between the accuracy of the method for obtaining airplane dynamic response as devel- 
oped in this study, and the accuracy of the methods based on the constant roll-rate assump- 
tion. In that comparison, the value of the constant roll-rate used was taken as the average 
value of the roll-rate variation p(t). In this regard, two questions naturally arise: 

1 . Is the average value of the rol l-rate the best value that can be used in the constant 
roll-rate methods? 

2. How does the accuracy of the method developed in this contract compare with the 
accuracy of the constant roll-rate method using values other than the average value of roll- 
rate? 

These questions were considered for the case of the family of rolling maneuvers studied 
previously. For each rolling maneuver, characterized by P m , the response in aand |3 was 
compared with the responses obtained from a number of constant roll-rate maneuvers ranging 
between 0.8 rad/sec to 3.0 rad/sec, and to responses 02 and $2 obtained by the methods 
used in this study. The results are shown in Table II . Each of the exact response in Qf(or 3) 
was compared with the approximate responses resulting from constant roll-rates of 0.8, 1 .0, 

1 .2, 1 .4, 1 .6, 1.8, 2.0, 2.5, and 3.0 rad/sec, and with the approximate response 02 
(or ^ 2 )/ obtained by the methods of this study. This is illustrated in figures 23 and 24 for 
the maneuver represented by P m = -2. 
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Table II indicates for each maneuver the values of constant roll-rate for which the correspond- 
ing response seems to give the least average error and the closest peak estimate. These are 
shown in rows 3 and 4 of the tabulation. Row 5 is the average of the values in rows 3 and 4, 
and may be considered as the "optimum" choice of constant roll-rate. Row 6 compares this 
optimum values of constant roll-rate with the average value ? QV of the roll-rate for each 
rolling maneuver. The following conclusions are evident from the tabulation and the figures: 

1 . The "optimum" value of constant roll-rate is 20 to 30% larger than the average 
value of the roll-rate for the particular family of maneuvers considered. 

2. Even if one uses this optimum value, the accuracy of the method of this work is 
far better than the accuracy of the constant roll-rate method. 

Analytical expressions for the time history of a maneuver starting from arbitrary 
non-zero initial conditions. - The method was extended to the case of a rolling maneuver 
starting from non-zero initial conditions in the variables Aa, |3, q, and r. The airplane 
response x(t) = [ Ac^(t) , p(t), q(t), r(t)] to an arbitrary rolling maneuver p(t), with initial 
conditions x Q = x(0) = [Aa^, q Q , r Q ], is governed by the equation 

x = A(p(t))x + u (p(t)) 

= A(P + p(t))x + u (p(t)) 

= A^x + A(t)x + u (p(t)) 


Repeating the steps of the method and taking into consideration the non-zero value of the 
initial conditions x(0), the successive approximations x n (t) can be expressed as 


t n 




W, (t,z)dz 


(33) 


where the terms W^(t,z) are defined by the following recurrence relationship: 


fH(t - z) [u(z) + 6(z)x ft ] 




\H(t - z) A(z) J W^_ 1 (z,y)dy 

0 


k= 1 


k = 2,3, . . .n 


A t 
o 


where H(t) = e and 6(.) is the unit-impulse function (delta function). 


(34) 
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TABLE II 


a 1 

3 

P m 

(Maneuver) 

1.0 1.5 2.0 

2.5 

i 

! 

3.0 

3.5 4.0 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

Average Value 
of Roll-rate 

.68 .87 .92 

i 

.94 1 

.945 

.95 . 955 

.68 

.87 

.92 

.94 

.945 

.95 

.955 

Pav 

1 











Value of Constant 
Roll -rate for 
least average 
error estimate 

■ 

.8 1.0 '1.0 

; 

1.0 ! 

1.0 

1.0 1 .0 

.8 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

Value of Constant 
Roll -rate for 
best peak 
estimate 

.8 h.2 1.3 

: } 

j 

1.3 | 

1.3 

* 

1.3 1.3 

.8 

1.0 

1.2 

1.2 

1.2 

1.2 

1.2 

Value of Constant 
Roll -rate for 
best "overall" 
estimate 

P c 

.8 1.1 1.2 

j : j 

| 

I 

j 

1.2 

! 

i 

1.2 

i 

1.2 1.2 

i i 

.8 

1.1 

1.2 

1.2 

1.2 

1 

1.2 

1.2 

P c 

Pav 

1.18 1.26 1.30 

1.28 

1 

t 

1.27 

1.26 1 .26 ’ 

! i 

1.18 

1.26 1.30 

i 

1.28 

; 

1.27 

1.26 

1.26 


i 




a 


c 

a \ 



FIGURE 3. COMPARISON OF EXACT AND APPROXIMATE RESPONSES 
IN a FOR ROLLING MANEUVER P m * -1 .5 RAD/SEC 
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FIGURE 8. COMPARISON OF EXACT AND APPROXIMATE RESPONSES 
IN a FOR ROLLING MANEUVER P » -4.0 RAD/SEC 




















-P (RAD/SEC) 
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FIGURE 16. COMPARISON OF EXACT AND APPROXIMATE 
PEAK VALUES OF a RESPONSE FOR FAMILY 
OF ROLLING MANEUVERS |P(t)j- p 





2.8 



FIGURE 17. COMPARISON OF EXACT AND APPROXIMATE 
FIRST PEAK VALUES OF p RESPONSE FOR 
FAMILY OF ROLLING MANEUVERS 


36 













FIGURE 23. COMPARISON OF a AND ou RESPONSES WITH 

RESPONSES RESULTING FROM VARIOUS VALUES 
OF CONSTANT ROLL-RATE (0.8 - 3.0 RAD/SEC) 




FIGURE 24. COMPARISON OF 3 AND ^ RESPONSES WITH 

RESPONSES FROM VARIOUS VALUES OF CONSTANT 
ROLL-RATE (0.8 - 2.5 RAD/SEC) 





Accuracy of the method for maneuvers involving prolonged durations of unstable 
ro 1 1 -rate. - The accuracy of the method was tested for rolling maneuvers that assume un- 
stable values of roll-rates for a relatively long period of time. The responses x^(t), X£(t) 

and x^(t) were obtained for the maneuver p(t) given by 


CN 

• 

T 

a> 

to 

• 

CN 

1 

II 

-4- 

^Q- 

0 < t <2.5 

(35) 

p(t) = -2.34e" 1,2t 

t > 2.5 

(36) 


Compared to the family of maneuvers studied in a previous section, this maneuver 
involves a relatively lengthy interval of unstable roll-rates, as shown in figure 25. Figures 
26 and 27 give the time histories for (3 and a respectively. As might be expected, these 
plots indicate that the approximate responses a n and |3 n exhibit a slower rate of convergence 
at the values of t for which p(t) is within, or close to, the unstable range of roll-rates. 
Nontheless, the agreement of the first approximations aj and with exact responses is 
substantially better than the agreement provided by a c and 8 C which are computed on the 
basis of the constant roll-rate approximation of previous studies. Table I II g ives the peak 
values of the various approximate responses: 


TABLE III 


a 

a 

c 


a 2 

a 3 

a 

a 

Oi 

a . 

a 

o 

o 

o 

o 

o 

1.02 

0.38 

0.95 

1.04 

1.12 


S c 



P 3 


P„ 

B 

p 

P 

o 

o 

o 

o 

o 

2.46 

0.16 

1.35 

1.92 

2.80 


Improving accuracy of method by piece-wise consideration of maneuver. - 
The convergence and accuracy of the responses * n (t) can be ameliorated by considering a 
given maneuver p(t) as a sequence of two maneuvers of shorter duration, with the terminal 
values of the first maneuver acting as the initial values of the second maneuver. This is 
feasible since we can obtain expressions for the time history of a maneuver starting from any 
initial values in the variables Aar, |3, q, and r, as given by (33) and (34). This procedure 
was applied to the maneuver considered in the previous section as given by (35) and (36). 

The responses x n (t) were first obtained for the initiation phase of the maneuver given by 

(35) , and the values x n (2. 5) were used as initial values for the termination phase given by 

(36) . Figures (28) and (29) show the time histories th us obtained for |3and a, which indicate 
a substantial enhancement in the accuracy (convergence) of the approximate responses over 

the corresponding responses given in figures (26) and (27), respectively. 
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Utilization of the Method - Usefulness of the Analytical Expressions 


In this section we illustrate the utilization of the method and demonstrate the usefulness 
of the analytical expression in the predictive study of the airplane dynamic response. Recall 
that the method provides analytical expressions for the mapping x(t) = M[p(t)], in which the 
response x(t) is expressed explicitly in terms of the maneuver roll-rate p(t). The analytical 
expressions for M may be useful in different ways. For any specific maneuver, one may sub- 
stitute the given expression of p(t) into M to obtain the response x(t), from which the peak 
excursions of the various variables may be obtained. The analytical expressions for the map- 
ping M may also be used to study the dependence (and sensitivity) of the response x(t) on the 
details of the maneuver p(t), such as level, rise time, slope, duration, etc. Simple ways for 
parametrizing the maneuver p(t) may be employed, and the mapping x = M(p) may then be 
used to obtain expressions for the response x(t) in terms of the selected parameters. Further- 
more, in view of the simpler analytical nature of the mapping x = M(p) in comparison to 
the original equations of motion, one may find it advantageous to use the expressions of the 
mapping x = M(p) to obtain the response x(t) on the digital computer with considerable savings 
in programming effort and computation time. 

To illustrate these points, a demonstration case-study was performed using the specific 
parameters of the F-100A. Explicit analytical expressions were obtained for the impulse 
response matrix H(t), which was then used to generate the approximate analytical expres- 
sions for a(t) and p(t) in response to the family of parameterized rolling-maneuvers (p(t)jp m 
described earlier. The analytical expressions were then employed to obtain information 
about peak excursions. Analytical expressions were obtained relating the magnitude of the 
response peak excursion to the value of the maneuvers index P m . All of the predictions 
based on the analytical expressions were in agreement with the information previously ob- 
tained via numerical integration. 

A t 

The impulse-response matrix H(t) .- The matrix H(t) =e appearing in the analytical 
expressions as given by (24) and (25) may be determined by any of the standard methods for 
obtaining the state-transition matrix of a linear system (ref. 9). The matrix H(t) was 
determined for the parameters of the F-100A as: 


H(t) = 10" 2 V(t) 


(71,11,29,8) 
(-100,84, .5,-66) 
(-23,-161,23,-81) 
(87,5, -87, -1.6) 


(3,47,3,-36) 

(29,10,71,-8) 

(105,8,-105,-8) 

(-15,50,15,99) 


(4,32,-4,11) 
(24,2, -24, -3) 
(67,-8,33,11) 
(-4,44,4,-55) 


(27,2,-27,-2) 
(68,-27,-68,-13) 
(.5,-17,-. 5,-50) 
(25,9,75,-12) 


where 


w/,\ _ / -»355t _ -,355t . _ -.206t q. -.206t , q,.\ 

V(t) - (e cos2.84t, e sin 2.84t, e cos.96t,e sin .96t) 

and where the multiplication V(a,b,c,d) is an inner product (dot product), e.g. 

V (71 ,11 ,29,8) = e" ,355t (71 cos2.84t + ll sin 2.84t) + e~‘ 206t (29 cos .96t + 8 sin .96t) 
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Analytical expressions for the approximate responses. - The matrix H(t) can be used to 
generate analytical expressions for the approximate time histories in response to a rolling 
maneuver p(t). Ignoring the effects of the weight terms and the elevator and rudder inputs 
(basic case), the input vector becomes 

u(t) = [0,p(t),0,0] (37) 


Substituting for u(z) and H(t - z) in the expressions (23) - (25), the following expression for 
a.(t) was obtained: 

a (t) = 10~ Z J p(t - z) [(-3,47,3,-36) V(z)]dz 

0 

Ignoring the smaller terms, one has 

a..(t) =*10 2 J* p(t - z)[47e " 3 ^ z sin 2.84z - 36e ' 2 *“*^ z s in .96z]dz (38) 

0 

Substituting the expression for the parametrized family of maneuvers described previously, 
namely, 

p(t) = P (1 - e”^*^) , 0 < t <i— H (Initiation) 

P m ” ' P rJ 

m 


= P 


m 



-0.6t 

e 


2tt (Recovery) 


the following approximate expression for the a response was obtained in terms of the 
maneuver parameter P^: 


Initiation Phase 

tt | (t)~ p fn [.226e' ,206t sin .96t + .2(e" ,6t -1) ] , 0^t<^I 

m 


Recovery Phase (t > 2rr/P m ) 




.226e"‘ 206t j^sin.96t-e 1 ' 3/ ' Pm 'sin.96(t-^-|) + .2(1 - e 3 ’ 8 ^ ? vi )e~ M 


Prediction of dynamic response information from the analytical expressions. - Considering 
these expressions for the range of values of P m of the family of maneuvers studied earlier, 
namely -A < P m < -1, the following information about the response in a was predicted in a 
straightforward fashion: 
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1 . The motion in a is a damped oscillation with a period of approximately 6 sec. 


2. The first relative peak value of occurs during the initiation phase of the 
maneuver at tj = 1 .4 sec, and the value of the peak is 

<* 1 (t 1 ) = a 1 (1.4) = 0.053P m 

3. The second relative peak value occurs at approximately + (2) (! 96) = 

2tt 

ip — . + 1 .65, i.e. 1 .65 seconds after the removal of aileron. The magnitude of 

' m' 

this peak is 


a ] (t 2 ) = P m [.tee" 1 * 3// ' Pm ' cos j- + ,07e" 3 * 8// ' ^-. 24 ] 

For the range of P m considered, lo^t^l > |ar^(t^)| , therefore, the expression for 
o^tJ represents the variation of the absolute peak value as a function of P m . Hence, this 
analytical expression corresponds to the graph shown in figure (30) which was obtained 
previously (see figure 16) from several time-history runs on the digital computer. 

It should be noted that the above information predicted by inspection of the analytical 
expression is in good agreement with the time-history runs in figures (2 ) thru (15) obtained 
by numerical integration of the equations of motion on the digital computer. 


Relating the maneuver roll-rate time history p(t) to aileron input 6 q (t). - It has been 
argued so far that the analysis relating the airplane dynamic response to the kinematic des- 
cription of the maneuver p(t), rather than to the aileron control input 6 a (t), is useful in its 
own right, particularily in a preliminary design effort where the characteristics of the flight 
control system are not fully specified. In the second portion of this study the airplane dyna- 
mic response to an arbitrary aileron input 6 a (t) was investigated. However, it should be 
noted at this point that for any rolling maneuver specified spatially by the kinematic des- 
criptor p(t), one may use the roll equation of motion (5) to obtain the corresponding aileron 
deflection 6a(t). In other words, equation (5) provides a means of deducing the aileron 
input 6 a (t) that would be required to realize a specified roll-rate time-history p(t). Equation 
(5) can be written as 




Substituting the expressions for (3, q, and r in terms of p(t), namely. 



one obtains an expression for 6 Q (t) in terms of p(t): 



M(p) M (p) - C M (p) - 

q r -cp p 


b 

2V 


(W )+ V)’ 


(39) 


Recall that the above relationships are valid under the assumption 



0 . 


Another consideration, which has a bearing on the problem of relating p(t) and 6 a (t), is 
the question of how to specify the shape of maneuver roll-rate time-history p(t), and the re- 
lated question of how to parameterized p(t) in meaningful and representative way for any 
given airplane. The results to be presented the remaining portion of this study on the air- 
plane response to specified aileron input ^(t) are helpful in answering these questions. A 
simple way to answer these questions for any particular airplane configuration would be to 
obtain a number of time-histories of the p(t) response to what might be considered as repre- 
sentative aileron inputs 6 a (t) applied during typical rolling maneuvers of the specific air- 
plane under consideration. 
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AIRPLANE MOTION IN RESPONSE TO AN ARBITRARY AILERON ACTUATION 


In this portion of the study, the dynamic response of the airplane resulting from a 
general actuation of the aileron input 6 a (t) is considered. Two types of results are obtained. 
The first type pertains to the derivation of analytical expressions of the airplane motion in 
terms of the aileron input 6 a (t), analogous to the expressions in terms of p(t) obtained in the 
first portion of the study. The second category of results pertains to the stability of motion 
under a general aileron input 6 a (t). Bounds on the motion variables are derived. A 
stability criterion is also presented. 


Analytical Expressions for Airplane Motion in Terms of Aileron Input 5 Q (t) 

The method presented in the first portion of this study was extended to relate the motion 
variables during a rolling maneuver to the aileron input 6 a (t) . The extended method 
generates analytical expressions for a, P, q, r, p in terms of an arbitrary given S Q (t), 
analogous to the expressions of a, P, q, r in terms of p(t) obtained previously. 

Outline of method . - Let the expressions for the motion variables in terms of p(t) be 
denoted by the mapping x = M(p): 


p(t) 



a 

M 

3 


q 


r 


The fifth equation of motion (roll equation), namely, 

p - Cp = C g & + BP + Dr + Eqr 
may be regarded as a mapping from (6 q ,x) into p: 


6 a 


3 

M c 

r 

lq 

5 


p=M 5 (6 o ,x> 


The iterative procedure for obtaining successively more accurate expressions for x and p in 
terms of 6 a is described as follows: For the given 6 a obtain from M5 a first estimate on p 
denoted as pj . Use p] into M to obtain an estimate xj = M(p-|), then use xj in M5 to 
obtain a second estimate P2, and so on. The accuracy and rate of convergence of this 
iterative procedure is dependent on the accuracy of the first estimate p-j . Two methods may 
be used to obtain the first estimate pj from the roll equation in a straightforward manner: 
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1 . Consider p^ as resulting only from the aileron input 6^(t): 

P, - C P, = C 6 6 a<’> 

2. Ignore the term Eqr in the roll equation and replace 0 and r by their expressions 

in terms of p(t) as obtained previously. The resulting integro-differential equation 

is a mapping from 5 (t) into p(t): 
a 

P - Cp - B0(p) - Dr(p) = C 6 § a 

from which an expression for Pj in terms of 6^ can be obtained. 

Derivation of motion analytical expressions in terms of 6 a (t) .~ Recall the expressions for 
the airplane response x n (t) = [a^t), 0 n (t), %(t)/ r n (t)] in terms of the roll -rate variation 
p(t). Denote the dependence of x^(t) on p(t) by the functional M^, 

x = [ a n /P'V r n J = M - n (p)/M(p)3 = M (p) (40) 

n n n n n an pn qn rn n 

where M (p), M. (p), M (p), M (p) are the expressions for a , 0 , q , and r in terms 
an r 0n r qn r rn r r n n t n 

of p(t) . We shall use these expressions together with the roll equation, namely, 

p -Cp -B0 - Dr -Eqr = C 6 6 a (41) 

to derive analytical expressions for a, 0, p, q, r in terms of an arbitrary given 6 (t). 
Rewriting (41) as a 


p - Cp = + B0 + Dr + Eqr = + f (42) 

where 


f = B0 + Dr + Eqr (43) 

The function f will be treated as a forcing function in equation (42) . For the given 6 a , one 
can use equation (42) to obtain a first estimate p-j (t) for p(t), as outlined in the previous 
section. Let a^, 0 n j, q n j, r n j be the expressions for a n , 0 n , q n , r n corresponding to p^ 
as given by (40): 

[a , , 3 w q ,, r .] = [M (p ), M (p.), M (p.), M (p,)] 
nl nl in nl an r l 0n r l qn r l rn 1 

Let f n l be the first approximation to f obtained from (43) by replacing 0, q, and r by 0^ , 
q^^, and r^, respectively: 
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where 


■ BP „1 + Dr nl + E %l r „l = %„<»!> + DM rn Cp,> + ^(p } ) MJp,) 
f nl<Pl ) = L n (p l> 



Using f n ] as an approximation to f, one obtains from (42) the second approximation p_ in a 
straightforward fashion: 


p 2- Cp 2 = C 6 6 a + f nl (p l>=Va + W 

p 2 = [C 6 6 a + L n <P l )] * e ° = J" CC 6 6 a (x) + " X> d * 

Repeating this procedure for P 2 to obtain p^, and so on, one obtains the following recur- 
rence relationships for p : 
r r m 

P m = [C 6 6 a + L n (p m-1 )] * ^ m = 2 ' 3 ' 4 ' • ' : < 44 > 

and the corresponding expressions for a, |3, q, and r 

a __ = M (p ) , 3 = M q (p ) , q = M (p ) , M (p ) (45) 

nm an r m nm pn r m ^nm qn r m ' rn r m 

The relationships in (44) and (45), together with the initial expression for p. , are the 
required analytical relationships representing the approximate responses in a, 3, p, q, and r 
for an arbitrary given 6 a (t). The index m represents the number of iterations and, therefore, 
the degree of accuracy, in relating p(t) to the given 6(t), while the index n represents the 
degree of accuracy in relating a., 3, q, and r to p(t), as described in the results of the 
previous sections which indicate that fairly accurate results should be obtained by taking 
n = 1 or n = 2. 

First approximation of p(t) in terms of 6 a (t) .- The initial approximation p 1 can be 
related to 6 (t) via equation (42) by either of the following methods: 

1 . Consider p. as resulting only from the aileron input 6 (t): 

I a 

P 1 - Cp l =C 6 5 a (t) 

p l ( f ) = C 6 6 a (0 * e Ct = / C 6 6 q (x) e C ^ f ^ dx 
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2. Ignore the term Eqr in equation (42) and replace 3 and r by their first approximate 
expressions 3-| and r^ in terms of p as obtained previously: 


p ) -Cp ) -B B| (p ) )-Dr 1 (p | )=C 6 e a 
P, - Cp, - Bp, * h 22 - Dp, * h 42 = C 6 6 a 


(46) 


where h22 and h^2 are entries from the impulse response matrix H(t) = [hjj] . The 
linear equation in (46) can be solved for pj in terms of 6 Q . Using the Laplace 
transform method one obtains an impulse response function h p § that relates p] to 6 a 
in (46): 

P l (t) = C 6 6 a * h p6 = C 6 ^ 6 a (x) h p6 (t ' x)dx 
where h - is the inverse Laplace transform of 

p6 


H p6 (S) S - C - BH 22 (S) - DH 42 (S) 

Method 2 provides a better approximation of 6 . However, method 1 involves a simpler 
computation . 


Stability of Airplane Motion Resulting from an Arbitrary Aileron Input 6 a (t) 

In the last part of this study, a method for analyzing the stability of motion in response 
to a general aileron actuation 6 (t) is presented. The predictive method provides bounds on 
motion variables. The method is also used to derive a stability criterion for boundedness in 
the form of an algebraic inequality in terms of the aileron input and the airplane parameters. 

Bounds on the motion variables . - Consider the roll equation 

p - Cp - BP - Dr - Eqr = C 6 = & 

o a a 

Replacing 3, r, and q by their expressions in terms of p, one has 

p - Cp - Bp * h 22 - Dp * h 42 - E(p * h 32 )(p * h 42 ) = 

where h 22 , ( 132 , h42 are the appropriate entries from the impulse response matrix H(t) = 

[hj - 3 . This nonlinear integro-differential equation relating p to 6 a is modeled into a feed- 
back system with a single product nonlinearity as shown in Figure 31 . The linear system G 
with a single input and two outputs is simplified in Figure 32, and characterized by the two 
impulse-response functions h and h : 

r q 
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( 47 ) 


r = (b -f)*h = £ *h - f * h 

' a r a r r 

q = (5 _ f) * h = 6 * h - f * h (48) 

^ a q a q q 


f = Eqr 


The functions 6 a * h r and £ a * h q , which represent outputs of the linear system G in response 
to the aileron input £ a (t), will be denoted by 6 Qr and 6 aq , respectively: 


6 = 6 * h , 

ar a r 

From equation (47) and (48), one has 


6 = 6 * h 

aq a q 


(49) 


In 5S l 6 ar l + l f * V "l 6 ar' + 


|q| ^|6 Qq l+ If * h q | -|6 aq l+ 


j f(x)h (t - x)dx 


j f(x)h (t - x)dx 


(50) 

(51) 


Consider the function f (t) defined as follows: 



Note that f (t) is a non -decreasing function, and that its increasing portions are identical to 
| f(t)| . Similarly, define the functions r(t), q(t), ^(O/ 6 (t) as 

r(t) = sup |r(t)| 

[ 0 ,t] 

q(t) = sup |q(t)| 

[0,t] 


6 

ar 


(t) = sup 1 6 (t)| 

[o,t] ar 


6 

aq 


(t) = sup 1 6 (t) 1 

[o,t] aq 


From (50) and (52) one obtains 


|r(0| ^ 6 Qr (t) + f(t) A r (t) 

(52) 

|q(0l * § (t) + f(t) A (t) 

(53) 


where 


t t 

A (0 s J |h (y)|dy = j |h (t - x)|dx 

0 0 r 

A (0 = J* |h (y)|dy = J |h (t - x)|dx 
q 0 H 0 s 


Note that A r (t) and A^(t) are non -decreasing functions of t. Since 6 ar (t), 6 a q(t), f (t) are 
also non -decreasing functions of t, it follows that the right-hand-sides of inequality (52) and 
of (53) are non -decreasing functions of t. Consequently it follows from (52) and (53) that 

sup |r(t)| S6 _ (t) + F(t) A (t) 

[o,t] ar 

sup |q(t)| ^ 6 (t) + T(t) A (t) 

[0,t] aq q 


or 


r(t) ^ 6 ar (t) + A f (t) f(t) 


q(t) * 6 aq^^ + A q^ 


And since f(t) = Eq(t)r(t) , one has 


hence. 


f(t) £ |E| q(t) r(t) 


r(t) s 6 (t)+|E|r(t)q(t)A W 
ar r 
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(54) 


Denoting |E| A^t) by a f (t) and |E| A q (t) by a q (t), 

r"(t) * 6 Qr (t) + a r (t) F(t) q(t) 

q(t)^6 (t) + a (t)F(t)q(t) (55) 

F(t) ^ 0 , q(t) 2 ; o (56) 

From these inequality conditions we shall now deduce conditions that guarantee the 
boundedness of r(t) and q(t), and derive explicit bounds R(t) and Q(t) on r(t) and q(t), 
respectively. Combining the inequalities (54), (55), and (56), keeping in mind that all 
terms are positive, one obtains 

r*(t)q(t) ^ [6 (t) + a (t) F(t) q(t) D C ^ (t) + a (t) r(t) q(t)] 
ar r aq q 


r (t)q(t) ^ 0 
Letting x(t) = r (t)q(t) one has 

a (t) a (t)x 2 (t)+ [a (t) 6 (t) + a (t) 6 (t) - l]x(t) + 6 (t) 6 (t) ^ 0 

r q r aq q ar ar aq 


Letting 


a(t) = a f (t) a q (t) 

b(t) = a (t) 6 (t) + a (t) 6 (t) - 1 

r aq q ar 


c(t) = 6 ar (t) 


Mt) 

aq 


the conditions may be expressed as 

a(t) x 2 (t) + b(t) x(t) + c(t) ^ 0 

x(t) s 0 


x(t) s 0 


(57) 

(58) 

(59) 

(60) 
(61) 


These two inequalities must be satisfied for all t ^ 0. We can now determine conditions on 
a(t), b(t), and c(t) which would ensure that x(t) satisfying (60) and (61) be a bounded func- 
tion . We shall assume that x(0) = r(0) q(0) = 0, which would be true for a maneuver starting 



from zero initial conditions. The analysis can be extended to the general case x(0) ^ 0. We 
shall also assume that x(t) is a continuous function. This is justified by the fact that the 
physical quantities r(t) and q(t) cannot be discontinuous variables. 


The inequalities (60) and (61), which must be satisfied for all t, impose a restriction on 
the values that x(t) may take. For any given t, let Xj be the set of all values of x(t) that 
satisfy (60) and (61): 


X 


t 


x(t) 


ax + bx + c ^ 0 
x ^ 0 


Figure 33 shows the set Xj for various conditions on the functions a(t), b(t), and c(t). Recall 
from (57), (58), and (59) that a(t) s0, c(t) sQ, but b(t) can take negative values. From 
Figure 33a, it is easily seen that if b^ - 4ac ^ 0, then Xj. = [0,°°3 . Figure 33b represents 
the condition b^ - 4ac >0, b > 0, for which also Xj = [0,*] . In Figure 33c , the condition 
b^ - 4ac >0, b < 0, forces the set X into two disjoint sets on the positive axis 

Xj = [0,^(0] U [r 2 (t),»] , r 2 > r ] 


where r^(t) and ^(t) are the roots of the equation 

a(t)x^ + b(t)x + c(t) =0 

Recall that x(t) G Xj for all t. Recall also that x(0) = 0 and x(t) is continuous . This implies 
that x(t) starts in the set [0,rj(t)] and remains in [0,r](t)] for all t. Thus, under the 
conditions 


b^(t) - 4a(t) c(t) > 0 , b(t) <0 , for all t 


(62) 


the time history x(t) is bounded by the function r^(t) 

x(t , S r ]W= ^V^ 


Recalling that x(t) = r(t) q(t), and substituting the expressions for a(t), b(t), c(t) from (57) 
through (59) one has 


r q £ 


1 - a 6 
r aq 


a 6 
q ar 


[1 - a 6 - a 6 

r aq q ar 


2a a 

r q 


- 4a 6 a 6 
r aq q ar 


Substituting into (54) and (55), and letting 


6 (t) =a (t) 6 (t) 
r q ar 


6 (t) s a (t) 6 (t) 

q r aq 


(63) 
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one obtains 


|r(t)|s r(t) * R(t) ^ 


1 - 6 (t) + 6 (t) -\/[l - 6_(t) - 6 (t) j - 46(t)6 (t) 


_g r 


q r 


q l 


|q(t)| ^ q(t) 5 Q(0 [ a 


1 - S (t) + 6 (t) -\/[l - 6 (t) - 6 (t)3 - 46 (t)6 (t) 


q" r" g r 

2a_(t) 


( 64 ) 


These expressions are the required bounds R(t) and Q(t) on r(t) and q(t), respectively. For 
any given aileron input 6 a (t) / one determines 6 r (t) and 6q(t) from (63) and (49) and substitutes 
in the above expressions to determine the time-history bounds R(t) and Q(t). Recall that the 
expressions for the bounds in (64) are valid only if the condition in (62) is satisfied for all t: 

[1 -6 (t) -6 r (t)] 2 -46 q (t) 6 r (t)>0 

6 (t) + 6 f (t) - 1 < 0 (65) 

These conditions will be used in a following section to derive a stability criterion which 
defines the range of aileron inputs 6 a (t) for which boundedness (non -divergence) of the air- 
plane motion can be guaranteed . 

Bounds on a(t) / P(t), and p(t) .~ Having determined the bounds R(t) and Q(t) on r(t) and 
q(t) , corresponding bounds on p(t), a(t) , and B(t) can be determined by relating these time 
histories to r(t) and q(t). From figure 31 it is easily seen that p(t) is the output of a linear 
system whose input is 6 a (t) - Eq(t) r(t). This is shown explicitly in figure 34, where the 
appropriate linear system is denoted by Lp, and its impulse response by hp(t). The time 
history bound P(t) on p(t) may now be expressed in terms of 6 (t), Q(t), R(t), and h (t) 
as follows: Q P 

p(t) = (6 - Eqr) * h (66) 

U 

lp(0l ^ |S (t) - Eqr | * |h | 

|p(t)l * C|6 o (t)| + |E| Q(t) R(t) ] . | h p (t) | ^P(t) 

The function P(t) is a time-history bound on |p(t) | . Alternatively, one may determine 
upper and lower bounds on p(t) which we shall denote as Pj(t) and P 2 (t), respectively. 

P 2 W 5 PW ^ 
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From (66), one has 


p(t) = 6 * h - (Ear) * h 
r a p ^ p 

|p - 6 * h | = | E | | (qr) * h | ^ |E| (|q| |r|) * |h | 
a p p P 

|p(t) - 6 (t) * h (0| « |E|[Q(t)R(t)l* |h (t)| 
a p P 

6 (t) * h (t) - |E |[Q(t)R(t)]* |h (t)| * p(t) £ 6 (t) * h (t) + |E|[Q(t)R(t)]* |h (t) | 
a p P a P r 

Thus, the upper and lower bounds on p(t) are: 

P^t) = 6 Q (t) * h p (t) + |E| [Q(t) R(t)] * |h p (t)| 

PJt) = 6 (t) * h (t) - |E| CQ(t) R(t)] * |h (t)| 

Zap p 

Having determined bounds on p(t), the corresponding bounds on a(t) and 3(t) can now be 
obtained from the expressions of a(t) and 3(t) in terms of p(t): 

a(t) = p(t) * hj 2 (t) , P(t) = p(t) * h 22 (t) (67) 

|a(t)| ^P(t)* |h 12 (t)| = n(t) 

|p(t)| * P(t) * 1 h 22 (t) | 

Alternatively, one may use the upper and lower bounds on p(t), namely P](t) and P2(t), to 
determine upper and lower bounds 0| (t) and ^(t) on a(t), $](t) and $2(f) on PW* F° r this 
purpose, we separate h -j 2(0 into two functions and hj 2 (t), re Pf esen Hng its positive and 

negative values, respectively, i.e. 

h 1 2 (t) = h 1 2 (t) whenever h 12 (0 ^0 

h ^ 2 (t) = h 1 2 (t) whenever h ^ 2 (t) £ 0 

Note that hj 2 (t) = h^ 2 (t) + hj 2 (t) for all t. From (67), one has 

a(t) = p(t) * [h 1 2 (t) + h 12 (t)] = p(t) * h 1 2 (t) + p(t) * h ]2 (t) 

and since P 2 (t) ^ p(t) ^ P ^ (t) , one obtains 

P 2 (t) h 12 (t) + P^t) h 12 (t) ^ a(t) s P^t) * h 12 (t) + P 2 (t) * h 12 (t) 
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Thus 


fij(t) = Pj(t) * h 1 2 (t) + P 2 (f) * h 12 (t) 

n 2 (t) = P 2^ * ^12^ + P l^ * 111 2 (t) 


Similarly for 8(t): 


$ 2 (t) s p(t) £ ^(t) 

where 

^W-P^t)* h 22 (t) + P 2 (t) * h 22 (t) 

^ 2 (t) ^ p 2 ( f ) * h 2 2 (t) + P l (t) * !l22 (t) 

Stability criterion. - In the two previous sections we derived expressions for the bounds 
on the time-histories of r(t), q(t), p(t), a(t), and 8(t). Recall that the existence of these 
bounds and the boundedness of the airplane motion was guaranteed if condition (65) is 
satisfied for all t: 


[1 - 6 (t) - 6 r (t)] 2 - 46 q (t) 6 r (t) > 0 

1 - 5 q (t) - 6^(t) >0 for all t (68) 


It con be easily verified from the definitions that 6q(t) and 6 r (t) are non -decreasing functions 
of time. This implies that the left-hand sides of the inequalities in (68) are non -increasing 
functions of t. Define 


d -mi (t . d - lim 6 (t) 
q a J r r 

“ f -»co ~ f- -*a> 


It follows that the conditions in (68) are satisfied iff 

[1 -d -d ] 2 -4d d >0 
q r q r 

1 -d - d >0 
q r 

or 

(d - 1) 2 + (d - I) 2 - 2d d > 1 
' q r q r 

d + d < 1 

q r 


(69) 


(70) 


( 71 ) 
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Thus, if for a given aileron input 6 a (t), the parameters dq and d r satisfy the conditions in 
(71), the response of the airplane is nondivergent. The set of values (dq,d r ) which satisfy 
(71) is obtained graphically in figure 35 by plotting the boundary curves whose equations are 

(d - 1) 2 + (d - l) 2 - 2d d =1 
q r q r 

d + d = 1 (72) 

q r 

The first boundary curve is in fact a parabola inclined at 45° and tangent to the axis at (1,0) 
and (0,1). The shaded region in figure 35 represents all the points (dq,d f ) which satisfy (71) 
and for which boundedness is guaranteed. By combining the various definitions that led to 
the definitions of dq and d r , one can express dq and d r in terms of the aileron time-history 
input 6 a (t) and the system functions hq(t) and h r (t) of figure 32: 


d r = |E| 

00 

/ |h(»|dt 
0 q 

sup 

X 

X 

J 6 a (t) h r (x - t)dt 

(73) 

d q - IEI 

00 

J* Ih (t)|dt 1 

0 

L 

sup 

X 

X 

f 5 (t) h (x - t)dt 

o a q 

(74) 


Stability of motion resulting from a step aileron . - This criterion for boundedness, pre- 
sented in the previous section takes a simple form for the special case where 6 (t) is a step 
input. Let a 


6 (t) = constant = 6 
a a 


Then from (73) and (74) we obtain 

d = |E| A m |6 | s y |6 | 

r ' q r a' q 1 a' 

d = |El A m |6 | *Y |6 I 

q 1 1 r q 1 a’ r a 1 


where 


A *J IMOIdt 


0 


a ■ J ih mi 

0 


m = sup 
* x 


J M*)dt 


= sup 


S MOdt 

A * 


(75) 

(76) 
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In the (d^ , d f ) plane of figure 35, equations (75) and (76) are the parametric equations of a 
straight line through the origin of slope Yr/Vq. The intersection of this line with the 
boundary of the shaded region gives the largest value of ] 6 a | for which stability is guaranteed. 
This value is obtained from (72), (75), and (76) as: 


7T - 77\ 2 
J 


max 


Y - Y 

q r 


q r 


(77) 


If Y ~ Y f - Y, then l^al max “ l/^Y, which is also the limit of the expression in (77) as Y c 


tends to Y . 

r 


7 ] 
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